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Abstract 


2019 novel coronavirus (2019-nCoV) caused severe, large-scale acute respiratory 
disease outbreak in Wuhan, China. The 2019-nCoV has spread to other regions and 
countries around the world, which is seriously threatening human health. There is an 
urgent need to develop drugs for the prevention and treatment of 2019-nCoV. 2019- 
nCoV nonstructural protein 14 (NSP14) carrying RNA cap guanine N7- 
methyltransferase and 3'-5' exoribonuclease activities could be a potential drug target 
for intervention. NSP 14 of 2019-nCoV shared 98.7% similarity with the one (PDB ID: 
Snfy) of acute respiratory syndrome (SARS) Coronavirus. Then, the 2019-nCoV 
NSP14 structures were modelled by using SARS NSP14 (PDB ID: 5nfy) as template 
for virtual screening. Based on the docking score, 18 small molecule drugs were 
selected for further evaluation. The compounds, including Saquinavir, Hypericin, 
Baicalein and Bromocriptine, could bind the N-terminus and C-terminus of the 
homology model of the 2019-nCoV Nsp14, thus providing as a candidate drug against 
2019-nCoV for further study. 


1. Introduction 


In December 2019, a large scale, severe acute respiratory disease named as “2019 
novel coronavirus (2019-nCoV)” occurred in Wuhan, China, and has already spread to 
other regions of China and other countries around the world in the following one month, 
seriously threatening human health. There is an urgent demand to develop drugs for the 
prevention and treatment of 2019-nCoV. Coronavirus NSP plays an important role in 
the virus’ genome replication and transcription !?, and it is generally conserved as an 
important functional protein in the coronavirus family. Among the family, NSP14 
protein has both exonuclease and methyltransferase functions, which is important for 
replication and transcription of SARS and other coronavirus, thus providing attractive 
target for drug designs >. 

The amino acid sequence alignment revealed that the NSP14 of 2019-nCoV shared 
98.7% similarity with the one (PDB ID: 5nfy) of SARS (Figure 1). Thus, the 2019- 
nCoV NSP14 structures were modelled by using SARS NSP14 (PDB ID: 5nfy) as 
template. The N-terminus and C-terminus of 2019-nCoV NSP14 were designated as 
active sites for screening drugs. A total of 7496 drugs obtained in the ZINC database 
were subjected to the binding screening. Among them, 2100 drugs were approved by 
FDA, 4264 drugs were approved by other regulatory agencies besides FDA (world-not- 
FDA) and 1132 drugs are undergoing clinical trials but not yet approved 
(investigational-only). The docking was carried out using AutoDock Vinal.1.2. Ten top 
compounds showed the lowest negative vina score in a range of -8.6 to -9.7 kcal/mol 
were selected from the N-terminal domain of homology model (Table. 1), and eight top 
compounds with lowest negative vina score in a range of -8.7 to -9.7 kcal/mol were 
achieved from the C-terminal domain of homology model (Table. 2). More importantly, 
the compounds, including Saquinavir, Hypericin, Baicalein and Bromocriptine, not 
only could bind the N-terminus and C-terminus of the homology model (Figure 1. A, 
B), but also could bind the N- terminal and C-terminal active pockets of the 2019-nCoV 


Nsp14 (Figure 1. C, D). 


2. Materials and Methods 
2.1 Docking method 


The SARS NSP 14 amino acid sequence was downloaded from the PDB protein 
structure database (PDB ID: 5nfy). The 2019-nCoV amino acid sequence (Accession 
number: MN908947) was obtained from database of the National Center for 
Biotechnology Information (NCBI). The homology of above amino acid sequence was 
aligned using ClustalW. Homology model of the target protein was constructed and 
optimized by Modeller9.18 using crystal structure of SARS NSP14 (PDB ID: 5nfy) as 
template. A total of 100 independent structures were constructed, and the one with best 
DOPE score was selected for further energy minimization by Amber. 

The ligands were downloaded from the ZINC database (FDA, world-not-FDA, 


investigational-only, http://zinc.docking.org/substances/subsets/). The 2D structure of 


the compound was then converted into the corresponding 3D coordinates using the 
Babel server (http://openbabel.sf.net). Then the model was converted to pdbqt format 
by prepare _receptor4.py script with assigning atom type and partial charge. All 
rotatable bonds in the ligand were set as flexible for flexible docking. Vinal.1.2 was 
used for molecular docking. The docking boxes were selected at the N-terminal 
exonuclease domain (aa: 62-290) and C-termimal transmethylase domain (aa: 291-527) 


of 2019-nCoV NSP14 respectively. 


2.2 Binding free energy calculation 


Each simulation system was immersed in a cubic box of TIP3P water with 10 A 
distance from the solute. The Na* or Cl was applied to neutralize the system. General 
Amber force field (GAFF) 15 and Amber ff14SB force field were used to parameterize 
the ligand and protein respectively. 10,000 steps of minimization with constraints (10 
kcal/mol/A2) on heavy atoms of complex, including 5,000 steps of steepest descent 
minimization and 5,000 steps of conjugate gradient minimization, was used to optimize 
each system. Then each system was heated to 300 K within 0.2 ns followed by 0.1 ns 


equilibration in NPT ensemble. Finally, 5 ns MD simulation on each system at 300 K 


was performed. The minimization, heating and equilibrium are performed with sander 
program in Amber18. The 5 ns production run was performed with pmemd.cuda. Based 
on the 5 ns MD simulation trajectory, binding free energy (AG) was calculated with 
MM/GBSA_ method according to the following equation: © AGcai=AH- 
TAS=AEvawtAEeletAGgbt+AGnp-TAS, where AEele and Agvpw refer to electrostatic and 


van der Waals energy terms respectively. AGgb and AGnp refer to polar and non-polar 
solvation free energies respectively. Conformational entropy (TAS) was not calculated 
for saving time. Besides, the ligands were compared based on the same target, so it is 


reasonable to ignore the entropy. 


3. Results and Discussion 
3.1 Docking results of Saquinavir against 2019-nCoV NSP14 model 


Saquinavir, as the first FDA-approved HIV protease inhibitor, has been used in the 
treatment of patients with human immunodeficiency virus (HIV) infection since 1995 
é, Our docking results showed that five of the hydrogen bonds involving ASP-273, 
ASN-252, ASP-90, and LEU-253 were maintained upon the binding of Saquinavir and 
N terminus of 2019-nCoV NSP14 (Fig. 1A). Meanwhile hydrogen bonds involving 
ASN-386, GLN-313, GLY-333 and THR-428 maintained upon the binding of 
Saquinavir and C terminus of 2019-nCoV NSP14 (Fig. 1C). Saquinavir could bind to 
the N- and C-terminal active pockets of the 2019-nCoV NSP14 (Fig. 1B, D). The recent 
study from a drug-target interaction deep learning model showed that Saquinavir can 
bind to 2019-nCoV RNA-dependent RNA polymerase to inhibit the enzyme activity’. 
Our simulation results showed that Saquinavir can bind two active sites of NSP 14, thus 


Saquinavir could be as a candidate drug against 2019-nCoV for further research. 


3.2 Docking results of Hypericin against 2019-nCoV NSP14 model 


Hypericin as a main ingredient in traditional Chinese medicine- Hypericum 


perforatum L. (St. John's wort) has been demonstrated activity against RNA viruses in 


vitro by inhibiting viral replication *. The present docking results showe that three of 
the hydrogen bonds involving ASN-252, GLY-93, and HIS-268 are maintained upon 
the binding of Hypericin and N-terminus of 2019-nCoV NSP14 (Fig. 2A). The six 
hydrogen bonds involving ASN-306, ARG-310, ASN-422 and LY-336 are maintained 
upon the binding of Hypericin and C-terminus of it (Fig. 2C). Hypericin can bind to the 
N- and C-terminal active pockets of the 2019-nCoV Nsp14 (Fig. 2B, D). Hypericin has 
been proven to have inhibitory effects on human hepatitis C virus (HCV) and human 
immunodeficiency virus (HIV)’. Combined the present study, Hypericin may have 
potential antiviral effect against 2019-NcoV. The traditional Chinese medicine- 
Hypericum perforatum L as main composition of Shuanghuanglian oral liquid has been 
widely used for the treatment of viral influenza. However, Shuanghuanglian oral 
solution has been suggested for the treatment of 2019-nCoV, triggering a huge crisis of 
public trust in Chinese scientists. We suggested that anti-2019-nCoV effects of 
Hypericin should be detected in cell culture models of 2019-nCoV infection. This will 
help us have a good understanding whether it is a good method to use Chinese 


medicines for the treatment of 2019-nCoV or not. 


3.3 Docking results of Baicalein against 2019-nCoV NSP14 model 


Baicalein, a flavonoid compound isolated from the root of Scutellaria baicalensis 
Georgi (Huang Qin in Chinese), inhibit viral replication of parainfluenza, influenza A, 
hepatitis B, HIV-1, and SARS coronavirus !™!?, The present docking results showed 
that six of the hydrogen bonds involving ASN-266, ASP-273, GLY-93, GLU-92, and 
HIS-268 were maintained upon the binding of Hypericin and N-terminus of 2019-nCoV 
NSP 14 (Fig. 2A). Four hydrogen bonds involving ASN-386, ASP-331, and GLN-313 
are maintained upon the binding of Hypericin and C-terminus of it (Fig. 2C). Baicalein 
can also bind to the N- and C-terminal active pockets of the 2019-nCoV NSP14 (Fig. 
2B, D). The previous study showed that Baicalein as a novel chemical inhibitor could 
inhibit ATPase activity of NSP13 protein of SARS coronavirus !°. The present data 


suggests that Baicalin may bind to NSP14 protein to exert anti-2019-nCoV activity. 


Therefore, we suspect that the anti-2019-nCoV activity induced by the Baicalein could 


be valuable for further study. 


3.4 Docking results of Bromocriptine against 2019-nCoV NSP14 model 


Bromocriptine, a specific dopamine receptor agonist for the hypothalamus and 
pituitary, has inhibitory effect on replication the Dengue virus with low cytotoxicity 
(half maximal effective concentration, EC50=0.8-1.6 uM; and half maximal 
cytotoxicity concentration, CC50=53.6 uM) '4. Moreover, Bromocriptine inhibited 
Zika virus protease activities and exhibited synergistic effects with interferon-a2b 
against Zika virus replications’. It is interesting to find that Bromocriptine can bind to 
the N- and C-terminal active pockets of the 2019-nCoV NSP14 from our molecular 
docking results (Fig. 4B, D). The present results showed that three of the hydrogen 
bonds involving ASN-104, ASP-273 and GLN-145 were maintained upon the binding 
of hypericin and N terminus of 2019-nCoV NSP14 (Fig. 4A). There was one bond 
involving THR-428 maintained upon the binding of hypericin and C terminus of it (Fig. 
4C). 


3.5 The calculation of binding free energy 

Based on the 5 ns MD simulation trajectory, binding free energy (AG) was 
calculated by MM/GBSA method. The calculated binding free energies of Saquinavir, 
Hypericin, Baicalein and Bromocriptine for the N-terminus of the homology model 
were -37.2711+3.2160, -30.1746+3.1914, -23.8953+4.4800, -34.1350+4.3683 
kcal/mol, respectively (Table 3), while the calculated binding free energies were - 
60.2757+4.7708, -30.9955+2.9975, -46.3099+3.5689, -59.8104+3.5389 respectively, 
when binding to the C-terminus (Table 4). Taken together, the results demonstrated that 


Saquinavir had the strong binding free energy. 


4. Conclusion 


2019-nCoV NSP14, a bifunctional enzyme carrying RNA cap guanine N7- 


methyltransferase and 3'-5' exoribonuclease activities could be a potential drug target 
for intervention. 2019-nCoV NSP 14 shares 98.7% sequence similarity with the 
corresponding one in SARS. Thus, the homology models of 2019-nCoV NSP 14 was 
structured for virtual screening. Based on the docking score, 18 drugs were selected for 
further evaluations. Four drugs (Saquinavir, Hypericin, Baicalein and Bromocriptine) 
could bind to the N-terminal and C-terminal domains of 2019-nCoV NSP 14. Combined 
the anti-viral function of above four drugs reported by the published literatures, we 
suggest the anti-2019-nCoV effects of above four drugs should be evaluated in the 


cell culture models of 2019-nCoV infection. 
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Figure and Table legends 


À PHE-146 


Figure 1. The binding model of Saquinavir against 2019-nCoV NSP 14. (A) Interactions between 
Saquinavir (cyan) and associated residues (off-white) in the N-terminus of the homology model for 
2019-nCoV; (B) Binding models of Saquinavir (cyan) in the 2019-nCoV NSP14 protein N-terminus 
pocket (white surface); (C) Interactions between Saquinavir (cyan) and associated residues (off- 
white) in the C-terminus of the homology model for 2019-nCoV; (D) Binding models of Saquinavir 
(cyan) in the 2019-nCoV NSP14 protein C-terminus pocket (white surface). Numbers 


accompanying dashed yellow lines represents the interaction distance (A). 
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Figure 2. The binding model of Hypericin against 2019-nCoV NSP14. (A) Interactions between 
Hypericin (cyan) and associated residues (off-white) in the N-terminus of the homology model for 
2019-nCoV; (B) Binding models of Hypericin (cyan) in the 2019-nCoV NSP14 protein N-terminus 
pocket (white surface); (C) Interactions between Hypericin (cyan) and associated residues (off- 
white) in the C-terminus of the homology model for 2019-nCoV; (D) Binding models of Hypericin 
(cyan) in the 2019-nCoV NSP14 protein C-terminus pocket (white surface). Numbers 


accompanying dashed yellow lines represents the interaction distance (A). 
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Figure 3. The binding model of Baicalein against 2019-nCoV NSP14. (A) Interactions between 
Baicalein (cyan) and associated residues (off-white) in the N-terminus of the homology model for 
2019-nCoV; (B) Binding models of Baicalein (cyan) in the 2019-nCoV NSP14 protein N-terminus 
pocket (white surface); (C) Interactions between Baicalein (cyan) and associated residues (off-white) 
in the C-terminus of the homology model for 2019-nCoV; (D) Binding models of Baicalein (cyan) 
in the 2019-nCoV NSP14 protein C-terminus pocket (white surface). Numbers accompanying 


dashed yellow lines represents the interaction distance (A). 
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Figure 4. The binding model of Bromocriptine against 2019-nCoV NSP 14. (A) Interactions between 
Bromocriptine (cyan) and associated residues (off-white) in the N-terminus of the homology model 
for 2019-nCoV; (B) Binding models of Bromocriptine (cyan) in the 2019-nCoV NSP 14 protein N- 
terminus pocket (white surface); (C) Interactions between Bromocriptine (cyan) and associated 
residues (off-white) in the C-terminus of the homology model for 2019-nCoV; (D) Binding models 
of Bromocriptine (cyan) in the 2019-nCoV NSP14 protein C-terminus pocket (white surface). 


Numbers accompanying dashed yellow lines represents the interaction distance (A). 


Table 1 —Ten drugs selected from the N-terminal domain of homology model 


Table 1 -Ten drugs selected from the N-terminal domain of homology model 


Drug name ID Data Affinity (kcal/mol) 
Hypericin ZINC000003780340 Investigational-only -9.7 
Bromocriptine ZINC000053683151 FDA -9.4 
Tanespimycin ZINC000100014666 Investigational-only -9.1 
Idarubicin ZINC000003920266 FDA -9.1 
Emend ZINC000027428713 FDA -8.9 
Baicalein ZINC000034114798 World-not- FDA -8.8 
Saquinavir ZINC000029416466 FDA -8.7 
Delavirdine ZINC000018516586 FDA -8.7 
Silibinin ZINC000001530850 Investigational-only -8.6 
Golvatinib ZINC000043195317 Investigational-only -8.6 


Table 2 —Eight drugs selected from the C-terminal domain of homology model 


Table 2 —Eight drugs selected from the C-terminal domain of homology model 


Drug name ID Data Affinity (kcal/mol) 
Hypericin ZINC000003780340 Investigational-only -9.7 
Olysio ZINC000164760756 FDA -9.4 
Sovaprevir ZINC000085537149 Investigational-only -9.1 
Celsentri ZINC000003817234 FDA -9.1 
Saquinavir ZINC000003914596 FDA -8.9 
Maraviroc ZINC000101160855 World-not-FDA -8.8 
Baicalein ZINC000034114798 World-not- FDA -8.7 
Bromocriptine ZINC000053683151 FDA -8.7 


Table 3 —The calculated binding energies of ligand to the N-terminus of 2019-nCoV 


NSP 14 


Table 3 -The calculated binding energies of ligand to the N-terminus of 2019-nCoV NSP14 


Energy” Saquinavir Hypericin Baicalein Bromocriptine 
AE vaw -52.9602+2.9999 -36.4737+4.0922 -36.8721+3.4155 -45.8461+3.1764 
AEele -128.6886+2 1.2732 -78.5578 £10.7496 -47.8615+12.4900 -119.9028+8.3707 
AGgp 151.4060+21.5835 90.4227+7.8130 66.0255+11.0420 137.397+46.4239 
AGnp -7.0283+0.3128 -5.5659+0.2085 -5.1872+0.2001 -5.7839+40.2553 
AGecal -37.2711+43.2160 -30.1746+3.1914 -23.8953+4.4800 -34.1350+4.3683 


* AEvaw = van der Waals energy terms; AEcie = electrostatic energy; AGg» = polar solvation free 


energy; AGnp = nonpolar solvation free energy; AGcal = final estimated binding free energy 


calculated from the above terms (kCal/mol). 


Table 4 -The calculated binding energies of ligand to the C-terminus of 2019-nCoV 


NSP 14 


Table 4 -The calculated binding energies of ligand to the C-terminus of 2019-nCoV NSP14 


Energy" 


Saquinavir 


Hypericin 


Baicalein 


Bromocriptine 


AEvaw 
AE ete 
AGgp 
AGnp 


AGcal 


-70.4383+4.1035 


-38.8487+7.5603 


57.7780+£6.3018 


-8.7666+0.3476 


-60.2757+4.7708 


-45.729+2.4822 


-58.4555+£12.1238 


78.5444+10.3832 


-5.3546+0.2317 


-30.995542.9975 


-48.647343.5522 


-192.8463418.1708 


202.1598 £16.7035 


-6.9761 0.1614 


-46.3099+3.5689 


-61.4659+2.9431 


-62.3583+7.0875 


70.6739+5.6693 


-6.6602+0.2480 


-59.8104+3.5389 


